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Comment on non-Gaussianity in hybrid inflation 
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In the literature there have been incompatible estimates for the amount of non-Gaussianity in 
hybrid inflation. In this note we point out the sources for the discrepancies and show that the 
results for the amount of non-Gaussianity in hybrid inflation obtained by two different methods, 
namely, perturbing Einstein equation to second order and the separate universe approach, indeed 
are compatible. This provides confidence in the methods themselves and in the actual computation 
of non-Gaussianities. 



Introduction. Recently there has been considerable in- 
terest in the possible non-Gaussian component of the cos- 
mological perturbations. Non-Gaussian perturbations in 
hybrid inflation were first estimated using consistent sec- 
ond order perturbation theory 1 in , where the formal- 
ism of Acquaviva et al. [2| was used. The approach is 
to perturb the metric and matter sides of the Einstein 
equation, and use the resulting equations to obtain the 
dynamics of required quantities. These quantities are 
then used to find out the second order curvature pertur- 
bation during inflation. 

This method leads to a set of equations and, even after 
simplifying (although motivated and not too constrain- 
ing) assumptions, the final expression for the second or- 
der curvature is complicated, containing several nonlocal 
and time-integrated terms. An order of magnitude of the 
result was estimated already in the original study |]|. It 
was later re-estimated in [2j, where the result seemed to 
be much larger. 

A completely different method, the SN formalism or 
the separate universe approach (see e.g. 3), for comput- 
ing the second order curvature perturbation in hybrid in- 
flation was employed by Lyth and Rodriguez in 5] . The 
result disagrees with the earlier estimates and implies an 
insignificant contribution from the transverse field, er, in 
hybrid inflation. 

There is a discrepancy between the results obtained 
with the two different methods, thus raising doubts on 
the validity of the methods. However, we argue that in 
both pj and Q the time evolution of certain quantities 
is not properly taken into account. In this brief com- 
ment we present a re-estimate of the original expression 
for the second order curvature perturbation from the a 
field. We show that when all the time evolutions are 
correctly taken into account the order of magnitude es- 
timate does indeed agree with the result obtained with 
the separate universe approach. The second order curva- 
ture seems to be proportional to the slow roll parameters. 
Such a small curvature alone would make the resultant 



non-Gaussianity unobservable, but according to |fj] there 
is a further suppression of the nonlincarity parameter in 
this particular model due to the uncorrelated |3{ nature 
of the a non-Gaussianities. 

Original computation. The original estimation of the 
non-Gaussianity was obtained in [IJ by extending the for- 
malism of 2] for two scalar fields. The approach is to 
expand the perturbations of the metric and the matter, 
which consists of two scalar fields, up to second order. 
These perturbed quantities are then used to write Ein- 
stein equations to second order. The curvature pertur- 
bation, defined for one scalar field in the first order as 
2] TZ = ip + H^-, is written for two scalar fields and 
expanded to second order; here if) is metric perturbation, 
H = a/ a is the Hubble parameter, 8<p is the inflaton 
perturbation, and (p is the time derivative of the back- 
ground value of the inflaton field. The Einstein equations 
are then used to obtain the evolution of TZ to second or- 
der, which in turn is then used to estimate the amount 
of non-Gaussianity. 

The analytic calculation of the evolution equations 
in the case of two scalar fields becomes complicated. 
To alleviate these difficulties two simplifying assump- 
tions are made in Q], namely, it is assumed that the 
background value o~o — and that the potential does 
not have any terms linear in a. The latter assump- 
tion means that <jq — indeed is a local minimum. 
The assumptions are well motivated and not too con- 
straining, and they clearly apply to many other mod- 
els in addition to hybrid inflation, whose potential is 0] 

V = V 



A ^.4 



(J- 1 + 



In fact, at 



1 This is during inflation. Complete second order perturbation the- 
ory connecting inflationary perturbations to the observed CMB 
anisotropics has so far not been developed. 



least some number of e-folds after horizon exit of the rel- 
evant scales, the form of the potential can be taken to 

2 2 

be V = Vq + ^f-c 2 + ,J ^-^ 2 i where m a and m are the 
effective masses of a and tp, respectively. 

The two constraints cause the second field, a, to com- 
pletely decouple from the first order Einstein equations. 
Its behaviour in the first order is only governed by the 
Klein-Gordon equation. At the second order, the con- 
tribution of a becomes completely additive, i.e. in the 
evolution equations there are no </?<7-mixing terms. This 
enables one to make use of the result of for the infla- 
ton contribution in Therefore, what is new in is 
the contribution coming from the transverse field a. 
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The expression obtained in the original paper for 
the contribution of the a field perturbations to the sec- 
ond order curvature perturbation is rather complicated, 
containing several nonlocal terms and several time inte- 
grated terms. It reads 



1 
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+4A- 1 d l (5a 1 d l S<r 1 y - 2{6& 1 ) 2 + m^cn) 2 
+(e - f])6HA- 2 d i (d k d k Sa 1 d l S<T 1 y 
+(e - r ] )HA- 2 d t d l (d k Sa 1 d k Scr 1 y 
-3A- 2 d i (d k d k 5a 1 d i Sa 1 y 

--A- 2 d i d i (d k 5a 1 d k Sa 1 y dt - A- 1 d l (6& 1 d i 5a 1 ) 

+3A- 2 d t (d k d k Sa 1 d l 6a 1 y + ^A- 2 d i d i {d k 5a 1 d k 5a 1 y 
+3eHA- 2 d l (d k d k 6<j 1 d l 5a 1 ) 
+ e -^A- 2 d l d\d k Sa 1 d k Sa 1 )} , 



(1) 



where A 1 is the inverse Laplacian, Alp = (8ttGn) is 
the reduced Planck mass, 8o~i is the first order perturba- 
tion of the transverse scalar field a, and e = -^(tjf^) 2 

and rj = MpyQ^ are slow roll parameters; dot denotes 
derivative with respect to time. 

In [| the slow roll solution 8a\ oc e^ m2 " t ^ 3H was used 
to obtain an estimate for the time derivative |<5<7i| ~ 

2 

^jf-\8<Ti\. (Double time derivatives were estimated by 
d 2 /dt 2 ~ (m 2 /H) 2 , but actually one should use equa- 
tion of motion to get rid of them) . Since both fields, ip 
and a, are effectively massless, the relation \8ai\ ~ \8ipi\ 
for the first order perturbations was used to approximate 



(2) 



In the first order Einstein equations a is completely de- 
coupled, and the situation is essentially that of a single 
field inflation. Therefore, "RP-' stays constant outside 
horizon. 

Only order of magnitude estimate was pursued and 
cancelling spatial derivative operators were neglected, 
e.g. lA^diK^d'H^l ~ \K^\ 2 . For the estimation of 
the time integral the quantities H , e, 77, m a and 8o~\ in the 
integrand were taken to be constants. The original result 
in pj for the estimate of the second order perturbation 
due to a reads 
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(3) 



iu pv aS^ 



the entire coefficient is of 



Since n a _ _ Py drj2 _ ff2 
the order slow roll parameters. 

Re-estimate by Lyth and Rodriguez. Later Lyth and 
Rodriguez [j| made a re-estimation of Eq. Q by inserting 



initial conditions and writing the equation as a definite 
integral. They, however, used the same estimates, Eq. © 

(with T^ 1 ) ~ const) and \5&i\ ~ ^f-\Sai\, as the original 
study 0. Similarly, they also assumed H, m a , and e to 
be constants, ending up with 2 



ftC0( t )_-fc(a)( t . 
+ml(Sa 1 ) 2 - 2(8a 1 ) 2 
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dt 



AN^\n^\ 2 



(4) 



Since AiV is the number of e-folds, which can very well 
be ~ 60, this result implies a much larger er contribution 
to the second order curvature perturbation than in the 
original study. 

Separate universe approach. In addition to the cos- 
mological perturbation theory approach 0, there also 
exists a recent computation ||| of the second order cur- 
vature perturbation in hybrid inflation using the separate 
universe approach. 

The general idea of the separate universe approach (see 
e.g. Q for a concise description) is to consider each point 
in space as being surrounded by a homogeneous FRW 
universe. Each point then has its own expansion param- 
eter N, i.e local number of e-folds, independent of the 
value of the expansion parameter (or any quantity) in 
other points. This expansion parameter depends on the 
values of relevant quantities, such as unperturbed scalar 
fields, at that point. The complete, inhomogeneous, be- 
haviour of the universe is obtained when all the sepa- 
rately treated points are patched together. 

The curvature perturbation, £, is in p| defined by 



9l] =a 2 (t)e 2 ^ %J (t,x) , 



(5) 



where, within inflationary context, 7^ contains the ten- 
sor perturbation which we do not consider here, (see Q 
for more details). Up to second order in scalar field per- 
turbations (8<j>i = 5<f>i(t,x)) C is obtained from [j| 



at, x) = j2 N ti {t)Hi + \y. 



(6) 



whereAT^f^andiV^^. 

Adapting the notation V = Vq(1 + \r\ip 2 + \rj fl cr 2 ) for 
the potential, the curvature perturbation reads |5j 

c= fy_Vrfyy + Va e 2AN(r,-r l „)(foy (7) 
■qip 2 \T)tpJ 2 \rjipJ 

The fields ip and a are assumed to be massless and 
their perturbations, Sip and 5a, are assumed to have 
the same spectrum (Ip-) in 5]. Therefore, we can set 



2 They used a different definition for the second order curvature, 
but that is irrelevant here. 
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1 ^ 1 


-0 — 1 


1 r]<p r 


1 i]tp \ 



ICil 

i.e. the contribution of a to the second order curvature 
perturbation as 



C2, CT ~ Oft*) e 



Cil 



(8) 



Source for the discrepancy. The seeming discrepancy 
between the two different methods basically comes down 
to articles and , or more precisely, assuming e to be 
constant, and to the assumption in Eq. i.e. 



|<M*)I ~ \^x{t)\ , 



(9) 



which is not generally justified. Eq. only holds im- 
mediately after horizon exit (t — U), when the amplitude 
of the perturbation of any effectively massless field / is 
\5f\~H. 

Using the slow roll equations (outside horizon) we ob- 
tain 



Sai{t) 

6<P!(t) 

<A)(0 



M(tOe-^ AJV , 
5^(U)e^ AN , 
MU)e^ AN , 



(10) 



where AN = HAt is the number of e- folds since ti\ thus, 
we obtain S&i = —r/ a H5ai, Sipi = —r/^HSipi, and ipo = 
—r/^HipQ. Since e ~ (p\j H 2 M^, we can also readily write 



e(t) = e ie - 2 ^ AN , 



(11) 



where we have denoted = e(ti). 

Now, it is immediately clear that \1Z^ \ = \H5<pi/(po\ 
stays constant, but one also sees that AN e- folds after 
horizon exit 



\H*^-\ ~ |7^>| & ~ e Aiv^-^)|^(i)| _ (12) 

tfO Otpi 



For the estimation of the time integral in Eq. the 
important point is that one may not move 1/e and So~i 
into and out of the time integral. 

Re-estimate of non-Gaussianity. Now we present a re- 
estimate of Eq. using the time evolutions expressed 
in Eqs. I|10|) and 1|11[1 . First, we notice that 

6HA- 1 d l (5a 1 d l 5a 1 ) + 2A~ 1 d l {8a 1 d l 5a 1 y 
-(S&rf + mKSo-i) 2 

= 2A- 1 d t [(3HSa 1 + 8a x + m 2 a 8a l )d i 8ai] 

= , (13) 

since 3H6&i + <5<7i + m^Sai — outside horizon. Thus, 



Eq. can be written 3 
1 



+2H(e - rj)% - %}dt - A^d^Sa^Sa,) 

+A- 1 d i (5a l d i Sa 1 ) + H(e - 277)7,} , (14) 
where we have denoted 



+ ^A- 2 d i d i {d k 6cj 1 d k 5a 1 ) 



(15) 



The last step is due to time derivatives within the time 
integral (77 and H are assumed to be constants). 

Eq. ((13*1) is used to get rid of the first term in Eq. QJ, 
-stti J 6HA~ 1 di(5&id l 5ai)dt, which would give too 
large a contribution to the estimate. Note that the order 
of magnitude estimate does not take into account possi- 
ble cancellations and, therefore, provides only an upper 
limit. However, the cancellations can be treated explic- 
itly, as is done here. 

For estimation purposes we also adopt the potential 
used in 5], V = V (l + ±77^ + t^ct 2 ). The slow roll 
parameters 77, and 77, are constants, and we also set H 
to be constant. The time evolutions of 8o~i, Sep,, tpo, and 
e are given by Eqs. (fTUfl and (jTT|) . 

Since the order of magnitude estimate anyway gives 
an upper limit, and since e < £j for any time t > ti, 
we replace e with except in the factor j. We again 
neglect cancelling orders of spatial derivative operators 4 , 
and put 5 a 1 = —q a H5ai. The estimate for the second 
order curvature perturbation, Eq. l|14fl. thus becomes 



eHM 2 



vlH 2 \Sa!\ 2 + e lVrT H 2 \Sa 1 \ 2 
+77 <T iJ|<5ai| 2 + e^lMI 2 + 77i?|<5ai| 2 } 
-J^-jo^) J VaH 2 ^ 2 dt 
+0(ei,r),r) a ) ff|foi| 2 } . (16) 



3 We use the indefinite integral here, instead of definite one with 
initial conditions, since the initial second order curvature, or at 
least initial non-Gaussianity, is supposedly small enough to be 
safely neglected Q. 

4 Because of this, both terms in 70- are effectively the same. 
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We have now two terms to evaluate, namely 



H\Sai 



eHMp 1 11 £l M2 



2 _ \5<Ti (U)\ 2 2AN(y-y a ) 
I , n,r2 C 



H 



5a! (ti) 



2AN{ V -Va) 



e 2AAf( t? -r, CT )| : ^(l)|2 _ 



(17) 



and 



eHM 



X w J H 2 V a\Sai\ 2 dt 



eMp 



\6<Ji(t, 



)_\_ e -2r,„AN _ e 2AAT(77 - 77 CT )|^(l)|2 _ /jgN 



eMf, 

Therefore, our final estimate reads 

TV? ~ 0{e u 77, 77.) e 2AN( ^ | 2 . (19) 

Comparing to the estimate in Eq. (JSJ one sees that the 
two results are of the same form and of the same order. 
Instead of 0(ei, 77, 7y CT ) Eq. (JHJ) only has a factor 0(ri a ). 
Here we have, however, provided only an order of magni- 
tude estimate and there still is a possibility for cancella- 
tions which may have been overlooked. It may be worth 
pointing out that if the nonlocal terms are discarded in 
Eq. I|14|) and only the first term is estimated, we obtain 
an 0(77,7) coefficient only. 

Discussion. The maximum number of e-folds for the 
observable scales is ~ 60 and observational limits for the 
spectral index of the curvature perturbation require 77 1 < 
0.01. It is therefore unlikely that the exponential factor 
e 2AAf(r)-r)„) wou ^ provide any significant enhancement, 
and the overall factor of the order slow roll parameters 
gives the magnitude of the result. 

The quantity measured in CMB experiments is the 
nonlinearity parameter /nl, which is related to TZ^ 2 \ 
but according to [|| it is highly suppressed in this sce- 
nario. The projected sensitivity (using first order per- 
turbation theory) for an ideal experiment is no better 
than /nl ~ 1 even including polarisation Q • It was later 
realized by Vernizzi || that the the second order cur- 
vature perturbation defined in has an artificial time 



evolution TZ^ ~ (2e — rj)TZ^ . However, even taking 
this small effect into account, it seems safe to say that 
the non-Gaussianity produced in the hybrid scenario by 
the a field seems to be too small to be observed unless 
some key aspects of the scenario are changed. 

There are still conceptual and practical problems with 
computing non-Gaussianities. The second order theory is 
not yet well established and the connection between the- 
oretical calculations and CMB observations is far from 
complete. One problem is that the usual scalar-vector- 
tensor decomposition of the perturbations is inherently 
non-local 0. These non-localities do not appear in the 
first order, but in the second order equations there are 
terms like lA~ 1 (d l g& L g) and A~ 1 (gAg), where g repre- 
sents a generic perturbation. The physical interpretation 
of these terms is not clear. 

Lyth and Rodriguez apply the separate universe ap- 
proach in [f| to compute the second order perturbations, 
or non-Gaussianity, in various scenarios including hybrid 
inflation. The formalism they use does not produce non- 
local terms involving the inverse Laplacian A -1 , and they 
state that ||| "... such terms must cancel if correctly eval- 
uated." Needless to say, this is quite a strong claim. It 
is actually not certain whether the separate universe ap- 
proach is completely correct when expanded to second or- 
der in perturbations. Indeed, in £| it was stated that any 
nonlinear interaction introduces mode-mode couplings 
which undermine the separate universe picture. As an ex- 
ample, recent studies of non-Gaussianities in preheating 
|l(j] demonstrate that these mode-mode couplings seem 
to be important. 

Despite all the problems and ambiguities in the dif- 
ferent second order formalisms it is comforting that the 
two different approaches discussed here can now be seen 
to produce the same result for the second order curva- 
ture perturbation and, therefore, for the amount of non- 
Gaussianity in hybrid inflation. 
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